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We investigate in detail gravitational waves in an Schwarzschild-anti-de Sitter bulk spacetime 
surrounded by an Einstein static brane with generic matter content. Such a model provides a useful 
analogy to braneworld cosmology at various stages of its evolution, and generalizes our previous work 
gr-qc/0504023 on pure tension Einstein-static branes. We find that the behaviour of tensor-mode 
perturbations is completely dominated by quasi-normal modes, and we use a variety of numeric and 
analytic techniques to find the frequencies and lifetimes of these excitations. The parameter space 
governing the model yields a rich variety of resonant phenomena, which we thoroughly explore. We 
find that certain configurations can support a number of lightly damped 'quasi-bound states'. A 
zero-mode which reproduces 4-dimensional general relativity is recovered on infinitely large branes. 
We also examine the problem in the time domain using Green's function techniques in addition to 
direct numeric integration. We conclude by discussing how the quasi-normal resonances we find 
here can impact on braneworld cosmology. 



I. INTRODUCTION 

Resonant states play a pivotal role in understanding 
many physical systems from drums to molecules to black 
holes [jj] . They often take centre stage because they con- 
tain the key observable features that are characteristic 
of the physical system, irrespective of the initial condi- 
tions that excite them in the first place. Braneworlds, 
it appears, are no different 0, 0]- Even though a sys- 
tem may be infinite and therefore unable to support a 
countable collection of normal modes, there may exist 
discrete complex frequencies which are nevertheless pref- 
erentially excited by initial data - so-called 'quasi-normal 
modes' (QNMs). The real parts of these frequencies are 
typically dependent on the length scale of the system, 
while the imaginary part tells us the rate at which en- 
ergy is lost to infinity in each particular mode. A recent 
investigation 2] has shown how the one-brane Randall- 
Sundrum (RS) model exhibits exactly this behaviour: 
the continuous Kaluza-Klein (KK) spectrum of gravitons 
has preferential masses with very short lifetimes. 

Such resonant phenomena may have a significant effect 
in braneworld cosmology. 1 As is well known, the effec- 
tive 4-dimensional formalism governing perturbations on 
the brane is not closed; that is, there are dynamical de- 
grees of freedom on the brane whose evolution is deter- 
mined by non-local bulk effects. Any complete picture 
of braneworld cosmology must include these degrees of 
freedom, but this is notoriously difficult to implement in 
practice. Limiting the discussion to one-brane models, 
we note that some progress has been made in certain 
special cases H El H , O ft IToll , b y us ing various approxi- 
mation schemes Um£[ EI III El El 113 111 EI , and 
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direct numeric solution of the problem in 5 dimensions 

HJIUillil. 

An alternative line of attack arises from realizing that 
the behaviour of perturbations will be dominated by res- 
onances of the brane-bulk system, which are precisely the 
quasi-normal modes. Given knowledge of the nature of 
these resonances, one should be able to make reasonable 
predictions about the effect bulk gravity waves have on 
the brane without knowing the details of how they were 
generated. This is especially useful since the correct ini- 
tial conditions for bulk gravity waves in such models are 
not really known. 

There are two aspects to the quasi-normal mode prob- 
lem for cosmological branes: the nature of the bulk, and 
the motion of the brane. The latter is tricky because 
of the non-separability of the equations of motion (the 
boundary condition in particular). However, a static con- 
figuration is relatively straightforward, and is a useful ap- 
proximation to a cosmological model provided the brane 
motion is 'slow'. 

A Friedmann-Lemaitre-Robertson- Walker brane may 
exist in a 5-dimensional bulk provided that it is a piece of 
Schwarzschild-anti-de Sitter (S-AdSs) We shall focus 
on the case of a positive curvature S-AdSs , with an Ein- 
stein static brane around it, and choose the 'inside' part 
of the bulk (i.e., including the black hole), while enforcing 
Z2 symmetry across the brane. In a previous paper [j| we 
have investigated tensor perturbations when the brane is 
pure tension. In this situation we found that gravity is 
'delocalized': gravity is sucked of the brane by the pres- 
ence of the black hole, removing the GR-like zero-mode. 
Moreover, the presence of quasi-normal modes was inves- 
tigated and were shown to be quite heavily damped, but 
with a mass-gap between the fundamental mode and the 
higher overtones. 

The pure tension condition, while useful for compar- 
ison with RS models, forces the brane to be located at 
the photon sphere of the bulk, providing only a glimpse 
of the possible range of behaviour relevant to cosmology. 
In particular, if the brane is moving it will pass, in princi- 
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pie, from the white hole region through the horizon and 
out to infinity (or maybe re-collapse). Hence we shall 
consider all brane locations outside the horizon. 

In this paper we show that the gravity waves are in- 
deed dominated by quasi-normal modes, the character of 
which depend strongly on the brane location. Roughly 
speaking, the closer the brane to the black hole horizon 
the higher the imaginary parts of the QNMs: this corre- 
sponds to the stronger gravity of the black hole draining 
gravitons of the brane more effectively. On the other 
hand, for branes located far from the black hole, or on 
small scales on the brane (provided the black hole mass 
is not too large compared to the AdS length scale), the 
damping of the modes is so small that the gravitons are 
effectively confined to the region between the brane and 
the photon sphere of the black hole: in effect, the quasi- 
normal modes become (approximate) normal modes, and 
gravity waves behave akin to two-brane scenarios. Thus, 
we find what we are looking for: brane signals are effec- 
tively 'combed' in the frequency domain to be composed 
of a discrete spectrum of very lightly damped harmonic 
modes - sharp evenly spaced spikes. Somewhat regard- 
less of the initial data, these spikes are the main feature 
of the signal on the brane. 

The layout of the paper is as follows: In Section [H] 
we discuss the Einstein static brane in S-AdSs, and ten- 
sor perturbations therein. In Section IIIII we consider a 
qualitative wave-mechanics analysis of the problem, be- 
fore going on to calculate QNMs for these models using a 
series solution for the master wave equation. This tech- 
niques proves to be computationally intensive for brane 
locations far from the horizon, so in Section IIVI we use 
an alternative approach for finding QNMs with small 
imaginary parts - which we call quasi-bound states. We 
then undergo a detailed study of how these states behave 
in the three dimensional parameter space which governs 
these models. We then discuss the utility of our results in 
the time domain in Section^using both analytic Green's 
function and numerical methods. We find the existence 
of coherent states which would look like bouncing gravi- 
tons to a brane observer; we also consider Hawking radi- 
ation, and its distorted spectrum on the brane. Finally 
in Section IVII we summarise our results and consider the 
implications for cosmology. 



Here, Rq is related to the ADM mass of the black hole 
while t is related to the (negative) cosmological constant. 
It is convenient to rewrite / as 
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When the solution is written in this way, it is obvious 
that there is an event horizon at R = R^. For simplicity, 
we will use dimensionless (t, r) coordinates defined by the 
substitutions: 



R 
T 



r x R h , 
t x i? h , 
ds 2 5) x J? h . 



(4) 



In terms of these quantities, the bulk line element looks 
like 



ds 2 s) = -/ dt 2 + f- 1 dr 2 + r 2 dnj, 



t 2 

(r 2 + 7 2 + l)(r 2 -l) 



7 



(5a) 

7T- (5b) 



In this representation, the bulk geometry is completely 
characterized by the ratio of the AdS length scale to the 
horizon radius of the black hole 7 and the horizon is 
always at r = 1. We shall call solutions with 7 < 1 'big' 
black holes and solutions with 7 > 1 'small' black holes. 
In order to reintroduce dimensions to any quantity in the 
(t, r) coordinates, one has to multiply by the appropriate 
power of the horizon radius R^. 



B. Einstein-static branes and their matter content 

An Einstein-static brane is introduced by identifying 
some r = rb as the boundary £0 of the 5-manifold and 
discarding the portion with r > rt,. The brane geometry 
is then simply that of an Einstein static universe: 



II. THE EINSTEIN-STATIC BRANEWORLD 
A. Bulk geometry 

The bulk geometry of the brane universes considered 
in this paper are given by the Schwarzschild-AdSs line 
element, which is conventionally written as: 



ds 2 5) = -f dT 2 + r 1 dR 2 + R 2 dtt 2 , 

Rq R 
R 2 I 2 ' 



1 



(la) 
(lb) 



ds {i) = -dr 



r 2 dQl 



(6) 



where r = f{r\,)t is the cosmic time. The normal to So 
is selected to point away from the brane and is obtained 
by evaluating the following at r = tv 



-r 1/2 d a r 



(7) 



With this normal, one finds that the extrinsic curvature 
of S is 



K„ h = h c h V r n„ = - 



1 df 
2/V2 dr 



u a u b 



f l/2 



'lab, 



(8) 
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where evaluation at r = r\, is understood. Here, we have 
defined the objects 



u a = -f 1/2 d a t, 
Kb = 9ab -n a n b , 
Qab = Kb + u a u b ; 



(9a) 
(9b) 
(9c) 



where u a is the future-directed normal to constant t sur- 
faces, h ab is the induced 4-metric on the brane, and q ab 
is the spatial 3-metric. 

The matter content of the brane is determined as fol- 
lows: First, assume a brane stress energy tensor of the 
perfect fluid form 



b = PtotU a Ub + Ptotqab- 



(10) 



The fluid quantities ptot and ptot refer to the total brane 
density and pressure, including any contribution from the 
brane tension ex. Hence, these can be further decomposed 
as 



Ptot 



k 1 a, 



Ptot = Pn 



K (T, 



(11) 



where k is a dimensionless 5-dimensional gravity-matter 
coupling constant, while p m and p m refer to the den- 
sity and pressure of 'ordinary' brane matter. Imposing 
I12 reflection symmetry across the brane yields (via the 
standard junction conditions): 

K S ab = -2(K ab -Kh ab ), (12) 

which is then equivalent to 

Kab = -K (jPtot + ^Ptot) U a U b - \nptotqab- (13) 

Comparison of (JSJ) and (|13fl yields the following consis- 
tency conditions for an Einstein static brane: 



1 df 



f 1 / 2 dr 

f X/2 



o / Ptot . Ptot \ 



K Ptot 



G 



(14a) 
(14b) 



It is possible to derive these conditions from the effective 
Friedman equation governing a general brane universe 
embedded in S-AdSs: 



a 2 36 



(15) 



where an overdot indicates differentiation with respect to 
cosmic time. A time derivative of (I15f) yields the Ray- 
chauduri equation for a, and setting a = a = yields (|14|> 
when the first law rf(ptota 3 ) = — Ptot^(a 3 ) is enforced and 
the identification a = is made. 

Equations l|14f> comprise a system of 2 equations in 4 
variables {7, rt,, Ptot>Ptot}- Hence, there is a 2-parameter 
family of Einstein-static brane solutions. For our pur- 
poses it is useful to regard the size of the bulk black hole 
7 and the brane position rb as free parameters, while the 



matter characteristics {ptotjPtot} are inferred from l|14fl . 
It is easy to show that for the restrictions 7 > and 
rb > 1, there is always a solution for {7, rb, Ptot>Ptot}- 
Furthermore, the effective equation of state of the brane 
matter is 



Wtot 



gtat = 3r b 4 - 7 2 - 1 + 2 7 2 r b 2 
Ptot 3(r b 2 + 7 2 + l)K 2 - 1) 



< 



(16) 



The inequality implies that a necessary condition for an 
Einstein-static brane universe is 2p to t + 3ptot < 0. In 
other words, the total brane matter distribution must 
violate the 5-dimcnsional version of the strong energy 
condition. This makes sense since the brane is essentially 
a static shell of matter around a black hole. The only 
way to maintain such a configuration is to construct the 
shell out of gravitationally repulsive material, hence the 
violation of the strong energy condition. Note that does 
not necessarily hold for the non-tension contribution to 
{ptot,Ptot}; i-e., it is possible to have {p m ,Pm} consistent 
with the energy conditions depending on the choice of a. 

Our principle concern is to understand how gravity 
waves change with brane location. If the brane were al- 
lowed to move, these conditions would no longer apply 
for a given brane position. 



C. Tensor perturbations 

We now consider tensor perturbations of the bulk using 
the general results of Kodama and Ishibashi [24[ special- 
ized to S-AdSs. Under such fluctuations: 



9ab 



g a b + Sq a b, u a Sq ab = n a Sq ab = 0, 



(17) 



where 5q ab is the perturbed spatial 3-metric. 

The SO (4) bulk symmetry allows a mode decompo- 
sition of 5q a b in terms of tensor harmonics on the unit 
3-sphere S 3 . Let i,j — 1,2,3 and Vi be the metric- 
compatible covariant derivative on 5* 3 . Then, tensor har- 
monics are defined by 



□T 



(ft) 



where □ = V*Vj, and 
'.2 



(ft) 



= TrT (fc) , 



k = L(L + 2) - 2, L = 1,2,3,. 



(18) 



(19) 



Now, if 9 l are suitable angular coordinates on S 3 , the har- 
monic decomposition of the metric perturbation reads: 



5q ab = d a 9 i d b 9 j J2 rl/2 <Pk(t,r)T, 



(ft) 

ij ' 



(20) 



When this is substituted into the linearized Einstein 
equations, on finds that each of the ipk satisfies a wave 
equation 



dt 2 



d 2 ^k 
dx 2 



-V k (r)ipk, x = x(r) 



du 

W) 



. (21) 
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FIG. 1: The potential in the master equation governing tensor 
perturbations in S-AdSs with no brane present and L = 1 

where x is the tortoise coordinate and the potential is 



also holds after perturbation. 2 Hence, 
5K a p = -\f 1/2 d r 5h af3 =>■ ^l% S 1ai3 = d r Sq a p. (25) 

Then, making use of l|14b(l to substitute for ftptot/ -1 ' 2 
and simplifying gives 

O = d r (r- 2 6q a0 ) r=rb . (26) 

It is worthwhile to contrast this with the boundary 
condition for perturbations 5h a p in a Randall-Sundrum 
model with scale factor A(y) and in the Randall-Sundrum 
gauge: d y (A~ 2 Sh al s) — on the brane. Finally, we use 
(12(it term-by-term in the expansion (|20|) to obtain the 
boundary condition satisfied by the master variable: 

= d r (r-^ 2 i: k ) , (27) 
which is of a mixed Neumann-Dirichlet type. 



V k (r) = f 



15 

4^2 



Ak 2 + 11 9(7 2 + 1) 



4r 2 



(22) 



The L — 1 potential is shown in Figure for several 
values of 7. We see that for large 7, there is a barrier 
near x = and an infinite wall at some finite x > 0. For 
small 7 the barrier is absent, but the wall is still there. 
The form of Vk highlights one of the curious properties 
of AdS space, namely r(x) — > 00 at a finite value of x. 
This means that, unlike asymptotically flat spacetimes, 
it is possible for gravity wave, electromagnetic, or other 
types of signals to propagate from finite r to r = 00 in 
a finite amount of coordinate time. The infinite walls of 
the potential curves in Figure ^ indicate the x position of 
spatial infinity, and in the absence of a brane one needs 
to specify boundary conditions there to have a well-posed 
Cauchy problem for the evolution of tpk- 

But we are interested in the case where there is a brane 
and the asymptotic region is not there, so we need the 
boundary condition on tensor perturbations at the brane 
location. To find this, we need to introduce 4-dimensional 
coordinates y a , where the a, (3... indices run over the 
directions tangent So. Then, 4-tcnsors are related to 5- 
tensors by 



'-a/3 



d a x a d x b T ab , etc. 



(23) 



In the background geometry, the fact that the radial off- 
diagonal metric elements g ra are zero imply 



III. QUASI-NORMAL GRAVITY WAVE 
RESONANCES 

We have seen in Sec. Ill CI that tensor perturbations 
of the Einstein static braneworld are governed by a one- 
dimensional wave equation Q21f) subject to a boundary 
condition l|27() at the brane. In this section, we solve this 
equation via a series solution in order to find the reso- 
nant modes of the potential, which are analogous the the 
quasi-normal modes familiar from black hole perturba- 
tion theory. 



A. Wave mechanics analysis 



If we substitute ipk(t,r) = e lu)t ^>ku{x) into (J2J, we 
obtain a time-independent Schrodingcr equation, 



(28) 



with energy parameter E 
boundary condition 



; 2 and a non-standard 



0=[r- s / 2 (x)*ku(x)]' r -_ 



(29) 



at the position of the brane. However, this boundary 
condition can be directly incorporated into the potential 
via an attractive delta- function. In Fig. |2 we sketch the 
resulting potential for several different brane positions in 
the case of a small bulk black hole. There are two qual- 
itatively distinct regimes characterized by the effective 
equation of state of the brane matter. If this matter has 



K a p = -y 1/2 d r h aP . (24) 

2 An alternative way of saying this is that the vector r a = is 
Since tensor fluctuations have n a 5g a b — 0, this formula parallel to n a both before and after perturbation. 



5 




FIG. 2: Tensor gravity wave potentials for various Einstein-static braneworlds with a small bulk black hole. We have taken 
7 = 10 and L = 1 for these examples 



a phantom equation of state wtot < — 1, we see that the 
brane is located in between the bulk black hole's event 
horizon and photon sphere. We call this the 'close brane' 
scenario. If the equation of state is more conventional 
— 1 < w t ot < — f , the brane is located between the pho- 
ton sphere and spatial infinity. Correspondingly, this is 
known as the 'far brane' scenario. In between the two 
possibilities is the pure tension case with iu to t = — 1 ; 
when the brane is coincident with the photon sphere. 

The qualitative form of the potentials in Fig.|21can give 
us valuable clues as to behaviour of gravity waves in the 
ES braneworld. We see that in all cases, the potential 
vanishes like e 2KX as x — > — oo; i.e., at the bulk black 
hole horizon. 3 Here, k is the surface gravity of the black 
hole: 



2dr 



7 



(30) 



As discussed in detail in Ref. |3j], this means we cannot 
find solutions for vf^ with lo 2 > that are localized near 
the brane. The vanishing of Vfe(l) means that the horizon 
has zero reflectivity; that is, it is a perfect absorber of 
gravity waves. One cannot construct a stable bound state 
under such circumstances because there will always be 
energy leakage into the bulk black hole. This means that 
there is no Randall- Sundrum like normalizable zero-mode 
for any Einstein-static braneworld — the bulk black hole 
essentially delocalizes the brane gravity. 

If not stable bound states, what other types of reso- 
nances could the potentials in Fig. |21 support? Well, the 
existence of the attractive delta-function suggests that 
we may be able to find metastable bound states. There 
are modes that are 'almost bound' to the brane, but are 
eroded in time by the tunnelling of the wave-function to- 
wards the bulk horizon. We can characterize such modes 



by a complex frequency with Imu > 0; i.e., they are ex- 
ponentially damped in time. Note that the larger the 
imaginary part, the higher the damping and the shorter 
the lifetime. In this sense, they are the equivalent to 
the familiar quasi-normal modes of black hole perturba- 
tion theory. Indeed, metastable bound states are defined 
entirely analogously to QNMs — namely, they are solu- 
tions of the wave equation satisfying 'purely outgoing' 
boundary conditions. In this context, this means that 
ipk satisfies both (|27l) and 



■j(t+x) 



as x 



(31) 



3 Though not shown in Fig. |2] the vanishing of the potential at 
the horizon is true for all values of 7. 



Because we are imposing two boundary conditions on 
a single ODE, it follows that the spectrum of QNM 
frequencies {uji} is discrete. Now, consider the situa- 
tion where a compact pulse of gravity waves strikes the 
brane and scatters. A well-known result it that the 
late-time behaviour of the scattered waveform will be 
dominated by contributions from QNMs (so-called quasi- 
normal ringing), independently of the initial data. It is 
for this reason that it is important to know the quasi- 
normal frequencies of the brane: they represent the gravi- 
tational resonances of the system and hence tend to dom- 
inate the gravity wave spectrum seen on the brane. 

As an important aside, we note that it is theoretically 
possible to find QNM solution with Jmu>i < 0. Such a so- 
lution will be spatially normalizable because ipk vanishes 
exponentially as x — > — 00. It will also be exponentially 
growing in time, and hence represents an instability of the 
system. Such an instability would not be unprecedented: 
It is well known that the 4-dimensional Einstein-static 
universe is unstable to homogeneous perturbations |25| . 
and that a Schwarzschild black hole enclosed in a finite 
cavity is unstable to polar perturbations 26] (assuming 
Dirichlct boundary conditions). This provides another 
motivation for finding the QNMs of this system; i.e., to 
determine whether or not it is stable against tensor per- 
turbations. 

Before we calculate the quasi-normal frequencies for 
the Einstein-static braneworld, we comment on the quali- 
tative features of the various small black hole cases shown 
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B. Series solution 

Only so much progress can be made by looking at the 
potentials in Figs. [5] and ED So we now turn our attention 
to an algorithm for actually calculation the quasi-normal 
frequencies of the ES braneworld. We begin by analyzing 
the master wave equation (|2f[) in the frequency domain, 
and transforming the radial coordinate to 



(32) 



FIG. 3: The tensor potential for a large black hole, far brane 
ES scenario with (7, L) = (0.5, 1). As in Fig.^J a = 3/(r b )/r b 
while loo = x(oo) 



in Fig. 13 The pure tension potential has been extensively 
studied in Ref. |3j], where it was found that the QNM 
spectrum featured frequencies whose real and imaginary 
parts were of the same order. How might we expect this 
to change for the close and far brane cases? For the far 
brane case, we see that the peak of the potential and the 
brane form a kind of finite potential well. It is plausible 
that gravity waves could be temporarily trapped in this 
well, but would eventually tunnel out to infinity. We ex- 
pect a set of QNMs to be associated with this effect, and 
their real parts should be roughly related to the width of 
the resonant cavity A. The trapping efficiency of the cav- 
ity will be directly related to the lifetime of the QNMs, 
and we will see in Sec. HVI that the efficiency can be very 
high indeed, resulting in quasi-normal mode frequencies 
with tiny imaginary parts. 

It is more difficult to guess the behaviour of the QNMs 
in the close brane limit. This should be considered as the 
very strong gravity regime, where the delocalizing effects 
of the black hole horizon are strongest. So we may expect 
that the imaginary parts of the QNMs to be much larger 
than the real parts. In this case, the potential is well- 
approximated by a simple exponential, which allows for 
an analytic solution that we discuss in Appendix lAl 

Finally, we comment on the large black hole case. As 
seen in Fig. ^ there is no clearly defined potential bar- 
rier for 7 < 2, so the far brane resonant cavity we dis- 
cussed above will not be present for large bulk black 
holes. The actual far brane potential with 7 small is 
shown in Fig.[3| Near the brane, the potential is propor- 
tional to l/(x — Xoo) 2 , where Xqo is the ^-coordinate of 
spatial infinity. The exponential fall-off of the potential 
as x — > —00 is unchanged from the small black hole case, 
and 14 smoothly interpolates between the two extremes. 
The lack of a resonant cavity implies that there will be no 
gravity wave trapping for far branes, but the qualitative 
features of the near brane and pure tension case should 
be the same as the small black hole scenario. 



This coordinate maps the horizon at r = 1 onto z = 
and r — 00 onto z = 1. Dropping the k index labelling 
the different tensor harmonics, we write 



(33) 



which converts l|21|l into an ordinary differential equation 
(ODE) of the standard form: 



(34) 



where P and Q implicitly depend on u>, 7 and L. We will 
attempt to construct a series solution of the equation, so 
we need to understand the pole structure of the contin- 
uations of P and Q into the complex plane. One finds 
that P(z) has simple poles at: 

z = 0, 2, and f ± i( 7 2 + f )~ 1/2 . (35) 

The poles of Q(z) are all of order 2 and are located at: 

z = 0, 1, 2, and 1 ± i( 7 2 + f )~ 1/2 . (36) 

From this information, we see that the horizon z = and 
spatial infinity z — 1 are regular singular points of the 
ODE about which we can construct series solutions. Our 
problem involves discarding the portion of the 5-manifold 
with r > rb, so we should attempt an expansion about 
z = 0. Our solution will have a radius of convergence 
greater than or equal to the distance to the closest pole 
to z — 0; i.e., the solution will be valid for z G [0, 1]. 

Following the standard procedure, we substitute the 
ansatz 



4>u(z) 



m=0 



(37) 



into the ODE. The coefficient of the lowest power of z 
yields an indicial equation with roots: 



e = ± 



2(2 + 7 2 )' 



(38) 



Switching back to the x coordinate, we see that the two 
possibilities correspond to 



e ±wx as x 



-00. 



(39) 
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Hence, the two signs correspond to waves travelling into 
or away from the black hole horizon. For a QNM solu- 
tion we need to choose the former, which is the same as 
demanding that our perturbation vanishes on the past 
black hole horizon. Once the sign of g is fixed, it is a 
straightforward but tedious exercise to obtain a ninth- 
order recurrence relation for the a m coefficients. 

We have yet to impose the boundary condition at the 
brane, which is 

4L(2b) ~ |n>^(^b) = 0, z h = l-l/r h . (40) 

To deal with this, we define an TVth order approximation 
to (j)^ via the partial sum: 

N 

f't^^V™ ( 41 ) 
m=0 

Substitution of this into the boundary condition yields 
an equation Vn = 0, where Vn is a polynomial of order 
2iV in iu>. Hence, only for a discrete set of 27V com- 
plex frequencies ui n will our approximate solution satisfy 
the brane boundary condition. Now, suppose that in 
the limit N — > oo the roots of the polynomial sequence 
{Vn}n=i converge to infinite set of discrete frequencies. 
Then, such a set represents a collection of resonant fre- 
quencies of tensor perturbations of the bulk-brane sys- 
tem. 

So, the problem of finding resonances for tensor pertur- 
bations of Einstein-static braneworld amounts to finding 
roots of a sequence of polynomials {Vn} that are sta- 
ble in the N — > oo limit. The practical implementation 
of this scheme is as follows: First, we fix all the input 
parameters of the problem rb, L, and 7. Then, we be- 
gin with a low order member of {Vn}, of degree 30 in ui 
say, and find all of its roots using Newton's method in 
the complex plane. Then, we use these roots as initial 
guesses for Newton's method applied to the polynomial 
of the next higher order. If Newton's method fails to 
converge after some number of iterations, we discard the 
root. But if it does converge, we use the answer as an ini- 
tial guess for the next order polynomial, and repeat the 
procedure. We end up with a finite number of roots as 
a function of N: u> n — u> n (N). The final step is to reject 
any of these roots that do not seem to be converging in 
the N — y 00 limit; i.e., roots for which the 'convergence 
factor' 

A n (7V)= \u n (N) - u n (N - 1)\, (42) 

does not become small for N large. At the end of the 
day, we are left with a set of resonant frequencies for our 
choice of rb, L, and 7. 

How big must we make N in order to have accurate 
answers for {o-i n }7 There is no entirely straightforward 
answer, but note that for far branes with ^ > 1, wc 
have Zb — ► 1. This means that we need to evaluate our 
series solution near the singular point z = 1 in order 
to enforce the boundary condition l|40|l . But we expect 



the series to be poorly convergent near a singular point, 
so in order to satisfy (|40l) accurately, N must be very 
large indeed. Hence, to find quasi-normal frequencies for 
7*b ^> 1, we need to retain more and more terms in our 
approximate answer <j>u ■ This is the principal practi- 
cal limitation of our method: the computational cost of 
computing {w n } increases dramatically with increasing 
rb- Indeed, in Sec. IIVI we will need to resort to other 
means to find the resonant modes for rb very large. 

C. Results 

In Fig. we plot the several smallest L = 1 quasi- 
normal frequencies of the brane for 7 = 20 and 7 = 
1/2. As mentioned above, the maximum value of rb wc 
can consider is limiting by computing power, so we limit 
ourselves to the range rb S (1,5). The roots of {Vn} have 
mirror symmetry about the imaginary axis, i.e., if w„ is 
a quasi-normal frequency, then — is also a solution. 
Hence, we only plot roots with Rew„ > 0. 

We highlight the key qualitative features: 

• In both cases, the QNMs are purely imaginary in 
the close brane regime, i.e., for rb — ► 1. These 
are overdamped, or evanescent modes , that show 
pure decay without oscillation. Also, they are 
evenly spaced in the limit, and the smallest root 
approaches zero. In Appendix^ we show analyti- 
cally that 

lim bj n = inn, n = 0, 1, 2, 3 . . . , (43) 

where k is the surface gravity (|30|l . As rb is in- 
creased, pairs of imaginary roots merge together 
and then move off on a trajectory with nonzero 
real part. 

• Also in both cases, there is a tendency for the fre- 
quencies to fan out and then migrate towards the 
real axis as rb is increased further, past the pure 
tension threshold. 

• The two cases exhibit quite different behaviour 
in the far brane regime. For a small black hole 
(7 = 20), the roots move towards the real axis 
very quickly, ending up with Reuj n much larger 
than Imw„. This means that the several QNMs 
are lightly damped in this limit. On the other 
hand, the large black hole (7 = 1/2) has roots that 
drift more slowly towards the real axis, resulting in 
moderately damped modes. The exception is the 
fundamental, defined as the QNM with the small- 
est frequency, which becomes lightly damped in the 



4 This is the common nomenclature in the theory of electromag- 
netic waveguides and other fields. 
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FIG. 4: QNM frequencies as a function of brane position; close brane roots are colored red, while far brane frequencies are blue 
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FIG. 5: QNM spectrum for close brane (rb = 1.05), pure 
tension (rb = 1-42), and far brane cases (rb = 10.00) when 
7 = 20 and L = 1 



far brane limit. Obviously, the reason for this dis- 
crepancy is that small black hole potential (Fig. [5J 
has a resonant cavity while the large black hole one 
(Fig. |2l does not. 

The extreme close and far brane behaviour will be dis- 
cussed in more detail in Sees. [X] an d If VI respectively. 

But first, let us consider the implications of the pat- 
tern of quasi-normal modes for actual gravity wave sig- 
nals on the brane. To do this, we fix 7 = 20 (i.e. a 
small black hole) and consider three separate brane po- 
sitions corresponding to extreme close, pure tension, and 
extreme far brane positions. The four smallest QNMs 
for each case are shown in Fig. |SJ We have plotted the 
imaginary parts of the modes on a logarithmic scale be- 
cause the damping varies by orders of magnitude from 
case to case. The right axis shows the half-life, defined 
as the time it takes for the modes amplitude to decrease 
by a factor of 1/2. The principal difference between the 
close and pure tension cases is that the Re ui n = in the 
former. The half-lives of the mode are similar, and less 
than about 10. However, the fundamental mode in the 
far brane case has tj/2 ~ 10 5 , which is extremely long 
compared to the other cases. In fact, all of the far brane 
modes shown have very small imaginary parts compared 
to the other cases. This types of behaviour are consistent 
with the speculation of Sec. 1111 Al — the resonant cavity is 
a efficient gravity wave trap. 

We now perform 'numeric scattering experiments' for 
each of the three cases shown in Fig. |SJ These consist 
of numeric solutions of the PDE (|21|l subject to specified 
initial data. In Fig. El we show the results of such in- 
tegrations at the position of the brane. The initial data 
for each of the three cases was identical Gaussian pulses 
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FIG. 6: Numeric brane gravity wave signals induced by a 
Gaussian pulse striking the brane when (7, L) — (20, 1) and 
rb = 1.05, 1.42, 10.00 from top to bottom. In all cases, the 
incident pulse has a width of 2 and is centered about xq = 
Xh — 20 at t = 0. The 'best-fit' values of u in the top two 
panels were found by matching a Re (Ae 1 "') template to the 
late-time signal, while the 'series' value comes from Sec. lUICl 
The inset for the last case shows the Fourier transform of ipk 
between t — 25 and 250 



incident on the brane at t = 0. The resulting signals 
reflect the characteristics of the QNM spectra seen in 
Fig. The close brane waveform is a purely decaying 
exponential with no oscillations. On the other hand, the 
pure tension signal takes the form of a sinusoid with a 
collapsing envelope. In both of these cases, the funda- 
mental mode's lifetime is relatively large compared to the 
first overtone, so we see signals dominated by one quasi- 
normal frequency. We can thus fit Ke(Ae' lu ' t ) templates 
to the signal to numerically determine the fundamental 
frequency luq. These are the 'best- fit' values shown in the 
upper two panels. Also shown are the frequencies deter- 
mined from the series solution shown in Fig. [SJ and the 
agreement between the two methods is excellent. 

The far brane signal in the lower panel is markedly 
different. We see in Fig. that three of the far brane 
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QNMs have lifetimes longer that the horizontal axis of 
Fig. El and the fourth's half-life is of the same order. So 
we expect contributions from all of the modes, which is 
exactly what is seen: a superposition of lightly damped 
QNMs. The signal is not adequately described by a sin- 
gle complex exponential, so we cannot find ujq by fit- 
ting. Instead, we plot the absolute square of the fast 
Fourier transform (FFT) of the signal between t = 25 
and 250 in the inset. The power spectra is peaked at the 
real parts of the four QNMs, showing that all of them 
contribute to the signal. What is misleading about the 
Fourier transform, however, is the nature of the peaks. 
We shall discuss this in detail in the next Section, but 
for now let us note the following: The peaks are located 
at the QNM frequencies; the signal near these goes like 
expi(w + iT/2)t, implying the power spectrum behaves 
like a Lorentzian, cx X/[(uj - zu) 2 + (r/2) 2 ]. Thus we 
expect the widths of the peaks in Fig. to be of order 
of the square of the imaginary part of the QNM: this is 
extremely narrow, and not easily resolvable numerically 
To summarize this section: We have found that when 
the brane is very close to the horizon, all the quasi-normal 
frequencies are purely imaginary. This makes the grav- 
ity wave signal a simple decaying exponential, as in the 
top panel of Fig. For intermediate brane radii, the 
real and imaginary parts of the QNMs have roughly the 
same magnitude, resulting in a 'classic' ringdown signal; 
i.e., the monochromatic damped oscillations seen in the 
middle panel of Fig. EI The far brane behaviour depends 
on the size of the bulk black hole. If it is large, the far 
brane behaviour is not much different than the pure ten- 
sion regime. But if the black hole is small, the resonant 
cavity between the brane and the photon sphere allows 
for a number of lightly damped QNMs. This results in 
a long-lived multi-frequency brane signal, as seen in the 
bottom panel of Fig. EJ 



IV. QUASI-BOUND STATES 

As mentioned above, the main drawback of our se- 
ries method of finding quasi-normal frequencies is that 
it becomes unfeasible as r^ — > oo. This is indeed un- 
fortunate, because as Fig. [5] demonstrates, some of the 
most interesting gravitational waveforms are associated 
with large brane radii. As seen in Figs. 0] and |31 this is 
due to QNMs with very small imaginary parts. We are 
naturally curious about the behaviour of these modes as 
rb becomes larger and larger, which motivates us to find 
a different method of locating the system's resonances 
when the imaginary part is small. Our techniques will be 
based on collision theory in ordinary quantum mechan- 
ics [53, and numeric solution of the wave equation in 
the frequency domain. Most interesting effects are asso- 
ciated with the resonant cavity between the brane and 
the bulk potential peak, so we will mostly highlight the 
small black hole case in this section. In Sec. IIV El we we 
indicate how the conclusions translate into the big black 
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FIG. 7: Numeric solution of the master wave equation in the 
frequency domain for a far brane, small black hole scenario. 
The plot also shows the definitions of various quantities ap- 
pearing in the main text 

hole scenario. 



A. The scattering matrix and the trapping 
coefficient 

In the beginning of Sec. IIII Al we wrote ipk = e wti &ku 
and hence transformed the master wave equation into a 
time-independent Schrodingcr equation l|28|l with energy 
parameter E = uj 2 and boundary condition (|29|l . Here, 
we consider the (numeric) solutions of this ODE with uj 
real. As in Sec. IIII Bl we will simplify our notation by 
dropping the wavenumber fc label from all quantities. 

In Fig. we show the results of a typical integration 
for a far ES brane, small black hole, and a given value 
of uj. As could have been deduced from the potential, 
the solution satisfying the brane boundary conditions 
takes the form of a simple sinusoid as x — » —00. We can 
write this asymptotic solution as 

* w « lA^e-^ie-*"* + S(oj)e iux ], (44) 

where S(uj) — e 2lS is the so-called scattering matrix and 
S = S(oj) is the scattering phase shift. 5 Presented in 
this way, ^ w is an explicit superposition of right and left 
moving waves. 

Now, we consider the continuation of S to complex 
values of ui. Suppose that the potential supports a QNM 
with Imw„ <C 1. Then, we require that S have a pole at 



5 In our case, the scattering 'matrix' is clearly a scalar quantity. 
This is because we are really looking at single-channel scattering, 
where there is no change of identity of the graviton as it scatters 
off the brane. Multi-channel scattering is common in nuclear 
physics, where the projectile and target can turn into different 
particles after the collision, which requires 5 to be matrix-valued. 
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u> = u> n in order to ensure that ^ is a purely outgoing 
solution. For simplicity, we assume that this is a simple 
pole. Also note that the scattering matrix is unitary 
SS* = 1 when to is real. Under these circumstances, the 
leading order Laurent expansion of 5* about its pole must 
be of the form El 



S(co) 



(45) 



where <5b g is the 'background' phase shift, which is an 
unspecified, slowly varying real function of u for our pur- 
poses. If we now write the nth QNM as 



TtlV. 



(46) 



we see that the total phase shift, with lu real and near to 
w n , is of the form 



5(u>) ~ <5bg(w) — arcsin 



n 



(47) 



Because T n is small, the second 'resonant' contribution 
implies that 5(lo) will vary rapidly about to — w n . In- 
deed, we expect that the phase will change by jt as w 
is varied across the resonance. The traditional way of 
identifying a resonance is to calculate 6(uj) for a given 
potential and look for sharp features in sin 2 S(tu). 6 Such 
lines will be centered about to — vu n and have widths 
proportional to r n , which means that long-lived QNMs 
give rise to narrow features. 

But is there a more physical quantity than sin 2 8{ui) 
to examine? Recall that QNMs are supposed to be 
metastable bound states, so we intuitively expect that 
^ to be localized near the brane for to w w n . We can 
quantify this expectation by defining the trapping coeffi- 
cient: 



( \ Ah 



(48) 



where is the magnitude of *$> u on the brane (see 
Fig. [7J. Any state with rj{uu) ^> 1 will correspond to 
the graviton being more localized near the brane than at 
infinity; this will be a 'quasi-bound' state. In fact, if our 
potential supported a stable bound state with Imw n = 0, 
we would have i](uj n ) — oo. 

Quantitatively, more can be said about rj{uj) near res- 
onance by writing 



Hf u (x)=Be[e iS ^\x)] ) 
where v]/^" 1 *' (x) is the purely outgoing solution 



vF 



>(*) 




(49) 



(50) 



Algorithms very similar to this (so-called 'resonance' methods) 
have be en used to identify lightly-damped QNMs of neutron stars 

MM- 



Here, R and 9 are assumed to be slowly- varying real func- 
tions of lu. This gives 

= r/ 2 H = i? 2 (w) sin 2 [<5(w) + 6{w) - tt/2]. (51) 



Now, suppose that the phases satisfy 

<5b g (w) + 9(lu) w 7r/2 or 37r/2, for 10 
then we have 



for 



(52) 



(53) 



Under this assumption, we expect £ res to be peaked about 
the real part of each lightly damped QNM. Furthermore, 
the shape of the peak is Lorentzian (or Breit-Wigner) 
with a half-width at half-maximum equal to the imagi- 
nary part of the mode. Of course, if 5b g + 9 is not ap- 
proximately 7r/2 or 37r/2 there is still a resonant feature, 
but with a more complicated lineshape. 

The only way to constrain these 'nuisance' phases is 
via direct calculation. In Fig. |H1 we plot £(w) for the 
(7,r b ,L) = (20,10,1) case. The spectrum features a 
number of discrete peaks, and the positions of the first 
few peaks match the real parts of the QNM frequencies 
calculated in Sec. 1111 CI for this case. We can make this 
more precise by subtracting a smooth baseline from £(lu) 
to remove some of the non-resonant part of the spectrum; 
i.e., contributions from the variations of i5b g and 9 with lu. 
The insets of Fig. |H1 show magnifications of the baseline- 
subtracted spectrum about the expected positions of the 
first three QNMs from the series solution; and the fre- 
quency range of each inset has been set to be 20 times 
the imaginary part. The shapes of the three lines are vir- 
tually identical and extremely well approximated by the 
Lorentzian template l|53|) . Hence, this example suggests 
that it is reasonable to take <5b g + 9 « 7r/2 or 37r/2 near 
an ES quasi-normal resonance. Furthermore, we see that 
the height of the peaks gives us roughly R(lu) 2 , after the 
baseline is subtracted. We have actually calculated the 
spectral function £(w) for many other cases, and found 
that the first few peaks in f always have a Lorentzian 
profile. 

So, we have a straightforward algorithm for estimat- 
ing lightly damped quasi-normal frequencies: First, find 
i{io) from the numeric solution of the ODE l)28[l. Then, 
subtract a baseline and find all of the features in £ res 
with a Lorentzian profile. The real and imaginary parts 
of the frequencies are simply related to the position and 
width of the those peaks. We will call resonances iden- 
tified in this fashion quasi-bound states (QBSs) for two 
reasons: First, since they maximize the trapping coeffi- 
cient they are, in some sense, brane-localized modes with 
Ab ^> Aoo. Second, they all have lmto n 1; i.e., they 
are extremely long-lived excitations. Of course, QBSs are 
merely a sub-class of QNMs, but we will see below that 
they dominate the phenomenology for far branes. 

Finally, we should stress that this trapping coefficient 
method really just estimates the QBS frequencies. This is 
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FIG. 8: The square of the trapping coefficient = rf(uj) 
in the (7, rb, V) = (20, 10, 1) case. The insets show the magni- 
fied views of the spectrum with a smooth baseline — shown in 
the main figure — subtracted. As described in the text, we 
expect the baseline-subtracted peaks to have a Lorentzian 
lineshape (1531 with w n and T n given by the series solu- 
tion of Sec. IIII CI We define the root-mean-square discrep- 
ancy between the normalized peaks shown above ft and 
the series prediction f s , assuming a perfect Lorentzian, as 
e = [/(ft - Zsfduj] 1 ' 2 / J £ B dLU. This leads to e = 0.011, 
0.0073, and 0.0072, respectively, which are very small dis- 
crepancies indeed 

because the procedure is prone to several uncertainties; 
including the non-uniqueness of our choice of baseline, 
deviations from Lorentzian profiles due to the running of 
<5b g , R and 9 with oj, the validity of the Laurent expansion 
(|45|l . etc. However, having said that, we find remarkable 
agreement between this 'trapping coefficient' method and 
the series solution in the (7, rb, L) — (20, 10, 1) case: 




0.247721 + 5.90 x 10~ 6 i (series) 
0.247720 + 5.88 x 10~ 6 i (trapping) ' 

0.6852 + 8.45 x 10~ 4 i (series) 
0.6853 + 8.51 x 10~ 4 i (trapping) ' 



Such accuracy is more than sufficient for our purposes, 
so we now proceed to use the 'trapping method' to cal- 
culate the QBS frequencies that the series solution has 
difficulties with. 




— > 

x 



FIG. 9: Spatial wavefunctions 9 n (x) of the first ff quasi- 
bound states in the (7, rb, L) — (20, 250, f ) case. 



B. Resonant cavity quasi-bound state 
wavefunctions 

Much intuition about the QBS solutions can be gained 
by examining their spatial profiles $™ (x) = 'i'^ (x) . We 
plot these spatial profiles in Fig. [5| for the (7, rb,L) = 
(20, 250, 1) case. The mode with the smallest frequencies 
is labelled with n — and is the fundamental, while the 
other modes are the 'overtones'. The logarithm of the 
potential V(x) is shown in the background. The vertical 
offset of each \I>™ corresponds to the energy parameter 
E n = zu^ defining the solution. We see that the overtone 
energies are greater than the potential within the reso- 
nant cavity, while the fundamental's energy only becomes 
greater than V in the asymptotic region. (Referring to 
Fig. this can be expressed as Wq < hi, while w\ > hi 
for n = 1, 2 . . .) As in ordinary quantum mechanics, the 
wavefunctions have concave curvature (\I/"/^/ < 0) when 
their energy is greater than the potential, and convex 
curvature (^f"/^f > 0) otherwise. Hence, the overtone 
profiles oscillate in the resonant cavity while the funda- 
mental shows monotonic decay. 

The integer n actually counts the number of nodes - 
places where the wavefunction crosses zero between the 
brane and the peak in V. This is also the number of local 
extrema — or 'loops' — in the wavefunction in this region. 
That is, 'J' 1 has one minima in the resonant cavity, 'J 2 has 
one minima and one maxima, etc. And the fundamental 
has none. This type of behaviour is strongly reminiscent 
of energy eigenfunctions in quantum mechanics, as in the 
radial eigenfunctions of the hydrogen atom, for example. 

The distribution of energies causes the n = mode 
to be the only wavefunction sharply localized near the 
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brane. Conversely, the overtone wavefunctions are dis- 
persed throughout the resonant cavity (and beyond for 
high-energy modes) . This type of behaviour generalizes 
to all of the other cases we have looked at. 

The key point here is that the 'quasi-bound state' 
moniker is well deserved for the lowest-lying modes, 
whose amplitude at infinity is vanishingly small. The 
inset shows a magnified view of the n = 1 wavefunction 
in the asymptotic region to emphasize that it is actually 
sinusoidal and nonzero. The shape of the overtone wave- 
functions suggests a tunnelling process where 'most' of 
the graviton is confined to the right of the barrier, but 
a small portion leaks through to infinity. This leakage is 
more pronounced for higher energy modes, corresponding 
to larger imaginary parts - so wider Lorentzian profiles - 
for these frequencies. 

Now that we have a fuller understanding of the trap- 
ping method and the quasi-bound states which it finds, 
we turn to a deeper investigation of the three dimensional 
parameter space of the model. We shall look at varying 
7"b, L, and 7 in turn. 
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C. Far brane limit: recovering the zero-mode and 
evenly-spaced overtones 

The trapping method allows us to further investigate 
the trends in the QNM frequencies shown in Fig. 0] to 
larger values of in the small black hole case (7 large) . 
In Fig. I1UI we show the real and imaginary parts of 
the quasi-normal frequencies when (7, L) — (20, 1) as a 
function of . The most striking feature of this plot is 
the disparate behaviour of the fundamental and overtone 
modes. While the imaginary parts of the overtones ap- 
proach constant values as rb — > 00, Imwo decays to zero 
as a power law. Hence, for infinitely large branes, we re- 
cover a single stable graviton mode confined to the brane. 

But is this infinite lifetime excitation the ES general- 
ization of the RS zero-mode? In 4 dimensions, tensor 
perturbations of an Einstein static universe would oscil- 
late with fre quen cy y 'L(L + 2)/Yb according to comov- 
ing observers [25j . Taking into account the redshift factor 
between the bulk coordinates and a comoving brane ob- 
server, this translates into a bulk frequency 



y7b£(£ + 2) 
WGR(7j^b) = - ■ (54) 

n> 

In Fig. 1111 we plot the fundamental QBS frequency for 
a range of (7, rb) parameters. Also shown as horizon- 
tal lines in the top panel is the GR frequency in the 
rb — > 00 limit. We see very clearly that for all 7, too 
approached ujqr asymptotically. Hence, not only do we 
have a mode with infinite lifetime in the large brane limit, 
the frequency of this mode matches what one might ex- 
pect from ordinary 4-dimensional theory. These are the 
essential characteristics of the standard zero-mode in the 



FIG. 10: Quasi-bound state frequencies as a function of brane 
position. Note that the fundamental frequency ljo tends to a 
real value as rb — > 00, while the imaginary parts of the over- 
tone frequencies are asymptotically constant and non-zero. 



RS scenario, so we conclude that as rb — > 00 we recover 
the zero-mode in the ES model. 

This is a very sensible result. We saw above that the 
n = wavefunction is distinguished from the others in 
several important ways, including being the only mode 
to be truly localized on the brane and nowhere else. So 
it is not surprising to find that its frequency shows a 
different dependence on tv Furthermore, in Ref. 0] we 
argued that the very existence of the bulk black hole was 
responsible for the delocalization of the zero-mode in the 
ES scenario. Here, we see that as the effects of the black 
hole at the brane position get smaller and smaller as rb is 
increased, we 're-localize' brane gravity. In other words, 
the gravitational field of the bulk black hole strips the 
zero-mode off the brane, but when the brane is far away 
this delocalization becomes more and more inefficient. 

Some of the features seen in Figs. E3 an d El can 
be understood using WKB-approximation techniques. 
In quantum mechanics, decay constants for metastable 
bound states can be approximated by the WKB tun- 
nelling amplitude across the potential barrier responsible 
confining the wavefunction. The classic example of this 
kind of calculation is Gamov's description of alpha-decay 
using elementary wave mechanics. As seen in Fig. |5J the 
energy of the fundamental mode lies beneath the reso- 
nant cavity; i.e., ujq < hi. Hence, we can view its decay 
as a tunnelling process from the delta- function at the 
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1281 to write: 



FIG. 11: Frequency of the fundamental QBS for varying 7 
and rb. The horizontal lines in the top panel indicate the GR 
prediction for the frequency of tensor perturbations in an ES 
universe as rb — » 00. Clearly, Wo approaches ojgr in the limit, 
suggesting that we recover a RS-like zero-mode for infinitely 
large branes. 



brane position. This gives rise to 



Imwo ~ exp 



; dx 



(55) 



Here, Xq is the 'classical' turning-point such that Wq = 
V(xo). Recall from Figs. and that the potential 
diverges at spatial infinity. The implies that as rb — > 
00, the rightmost portion of the integrand diverges like 
(a^oo — icb) 1 , which in turn gives Imwo — > 0. Further- 
more, it is not difficult to show that 



Imwo oc r. 



/15/2 



as rb — ■> 00. 



(56) 



This is exactly the asymptotic behaviour seen in Fig. 1111 
log 10 Im u>o versus log 10 rb is a straight line for large brane 
radii whose slope is independent of 7. We expect this 
behaviour to generalize to other choices of L. 

Turning our attention to the overtone modes, we note 
that their energies satisfy hi < w 2 ^ < h,2, so there are 
three turning points xq < x\ < x-i. The potential barrier 
that is responsible for the quasi-confinement is the peak 
in V(x) shown in Fig. |5J This peak does not become 
infinitely large in the limit, which explains why the over- 
tones have finite damping asymptotically. In this case, 
WKB techniques can tell us something about the real 
parts of the QBS frequencies. Because the overtone en- 
ergies are greater than the potential inside the resonant 
cavity, we can use the Bohr-Sommerfeld quantization rule 



dx y/ — V(x) = (n + |)7r, n > 1. 



(57) 



If we approximate the potential in the cavity by its av- 
erage value, we obtain: 



ZD. n 



<v 



(n + i)7T 



(n + ±)tt 
Ax 



for V < 



(Ax) 



(58) 



where Ax = 22 — %i ■ The last approximation holds rea- 
sonably well for the higher modes in the current prob- 
lem. Finally, as rb — > 00 we have Ax ~ d, where d is 
the distance between the brane and the potential peak 
(cf. Fig. UJ. In the limit d becomes constant, which ex- 
plains why Rew n becomes asymptotically constant and 
evenly spaced for the overtones in Fig. ITU1 



D. Small scale limit: the quasi-bound states 
become bound 

So far we have concentrated on perturbations with 
L = 1; now consider perturbations with L not necessarily 
equal to unity. The value of L controls the wavelength 
of the gravity waves tangent to the brane, so when we 
increase L we are examining smaller and smaller scales. 
What kind of modes will exist in this case? 

In Fig.^1 we plot the logarithm of the trapping coeffi- 
cient £(u>) for the (7, r b ) = (20, 250) case, but now letting 
L increase from 1 through 5 to 10. The principle effect of 
increasing L is to increase the number and sharpness of 
the spikes - implying more quasi-bound states with much 
lighter damping. The WKB approximation can again be 
used to understand this behaviour. 

Returning to (J5SJ, we take the L — > 00 limit instead of 
rb — > 00. The height of the peak in the potential scales 
like 



(7 2 + 2) 2 2 

4 7 2( 7 2 + 1) ' 



(59) 



so its contribution dominates the integral. The same 
holds if we apply the formula to the overtones by switch- 
ing vjq to w n and Xh to x\. In any case, this leads to 



Imw„ cx e 



- a L 



as L 



(60) 



where a is a constant depending on 7 and rb. This ex- 
ponential decay in the widths of the spikes is shown in 
Fig. where we plot the imaginary parts of the QBS 
frequencies for various scales. 

Physically, this is intriguing: As the spatial scale of 
perturbations on the brane decreases, their bulk angu- 
lar momentum around the black hole increases. This 
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FIG. 12: Variation of the trapping coefficient with L. 
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FIG. 13: The variation of the imaginary parts of QBS fre- 
quencies with L. 



causes the centrifugal barrier between the brane and the 
event horizon to grow taller and taller, which in turn 
makes it difficult for gravity waves to tunnel away from 
the brane into the black hole. The net result is that the 
confinement of the QBSs in the resonant cavity becomes 
100% efficient. Effectively, the quasi-normal modes be- 
come normal modes when their imaginary parts disap- 
pear; the quasi-bound states become bound when 
goes to zero. This is like a two-brane model where the 
potential barrier takes the role of the shadow brane. We 
will discuss the implications for cosmology below. 




FIG. 14: £(u>) for various L = 1 big black hole scenarios 



How many QBS states are there in the large L limit? 
Above we saw that high overtones have u) n ~ 7r(n+ \)/d. 
The highest QBS mode should have have an energy com- 
parable to the height of the potential peak; i.e., Lojq ~ hi- 
Finally, we know that the height of the peak scales as 
L(L + 2), which yields that the number of quasi-bound 
states is proportional to L as L — > oo. 



E. The big black hole case 

To this point, we have only applied the trapping coef- 
ficient method to small black hole cases featuring a res- 
onant cavity. But under certain circumstances, we can 
also use the procedure when 7 is small. 

The square of the trapping coefficient for range of rel- 
atively small values of 7 and r b are shown in Fig. 1141 
All the spectra show one peak indicating the existence of 
only one QBS solution. This is because there is no reso- 
nant cavity when 7 is small (see Fig. hence there are 
no overtones. The general trend is for the peak position 
to show little variation with r\,, and to decrease with in- 
creasing 7. On the other hand, the width of the feature 
decreases as either 7 or r b is increased. 

The rather broad nature of the peaks for small (7, r b ) 
means that they are not well approximated by Lorcntzian 
profiles, which suggests no modes with tiny imagi- 
nary parts. This is consistent with the series results 
of Sec. 1111 CI where small imaginary parts were only 
achieved for large brane radii. When r b is large, the 
fundamental-mode WKB analysis of Sec. IIV CI is appli- 
cable to this case, implying that the imaginary part of 
the sole QBS frequency goes to zero for large branes. We 
have confirmed this numerically, and also confirmed that 
the real parts approach u>gr- So, we also recover the GR 
zero-mode as r b — > 00 in the big black hole case. 

As in the small black hole case, the potential barrier 
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between the brane and the horizon scales as L(L + 2) 
for large L. Hence, the conversion of quasi-bound into 
bound states as L — > oo also occurs for small 7, the chief 
difference being that there is only the fundamental mode 
to covert. 

The main difference between the small and large black 
hole cases, then, is that in the former there are quasi- 
bound state overtones, while in the latter there are not. 
However, there is always a fundamental n = QBS, and 
the limiting behaviour of that mode as or L becomes 
large is the same irrespective of 7. 



the boundary condition at the brane, x < = min(a;,a;'), 
X> = max(x,x'), and W = ^ L) *L R) ' ~ ^S '*?' is thc 
Wronskian. In this section, it is useful to choose the in- 
coordinate such that Xb = 0. Using the definitions and 
results of Sec. W~ 



*L R) (*) 



we have 



Re 10 , . 
cos(wa; - 



> —00, 
x b = 0, 

x — > - 



> —00, 
x b = 0. 



(64a) 



(64b) 



V. THE INITIAL VALUE PROBLEM 

In Sec. IIII CI we solved the wave equation in the time 
domain to illustrate the effects of the QNMs on actual 
brane signals. Here, we examine this problem in more de- 
tail by considering several different classes of initial data. 
We will demonstrate that the trapping coefficient 77 is 
useful for more than just finding QBSs; it actually plays 
a prominent role in determining the power spectrum of 
the brane signal. We also consider the scattering of wave 
packets in scenarios featuring a resonant cavity, and ob- 
serve coherent 'bouncing' behaviour. Finally, we allow for 
'initial data' representing thc situation where the brane 
is in thermal equilibrium with the semi-classical Hawking 
flux from the black hole. 



A. Green's function analysis 

Suppose that at t — 0, we specify the both the value 
and first time derivative of ip. Then, we can construct 
a formal Green's function solution for the signal on the 
brane ipb(t) = t/j(t,rt,) by following the standard tech- 
nique from black hole perturbation theory p], [^J . This 
involves applying the Laplace transform to the master 
wave equation to convert the initial data into a source 
term, obtaining the Green's function solution, and then 
performing the reverse transform to return to the time 
domain. If we write the resulting solution as an integral 
in Fourier space, we obtain 

ip(t,x) = — [ duj dx'e^Gufax'yiix' ,u), t > 0. 
2n J 



Here. 



I(x', uS) = -iojip(0, x') - ip(0, x'). 



(61) 



(62) 



As usual, the Green's function is composed of solutions 
to the frequency domain homogeneous equation l|28(l : 



G LJ (x, x ) 



(63) 



where is the ODE solution that satisfies the out- 
going boundary condition at x — —00, ^JJ 1 ' satisfies 



Because there is no linear derivative in the Schrodinger 
equation l|28|l , the Wronskian is independent of x and we 
can evaluate it at x = — 00 



(65) 



If we now select the field point to lie on the brane r = r b , 
our solution for ip reads 



where 



1 

2^ 



du;e lut e iS r](u:)i(Lo), t > 0, (66) 



dx^ ] {x) 



4(o,x) 



(67) 



Note the prominent role played by both the scattering 



matrix S 1 ^ 2 



and the trapping coefficient rj in l|66|l . 



There are two ways one can evaluate the ui integral: ei- 
ther by direct or contour integration. In the former, 77 
appears directly in the Fourier transform of the brane 
signal. Since rj is sharply peaked near QBS resonances, 
the frequency profile of the brane signal will be similarly 
peaked. On the other hand, if one does the integral by 
completing the contour in the upper-half plane, the poles 
in the scattering matrix will give the dominant contribu- 
tions to the signal; and of those contributions, the QBS 
poles will be the longest-lived and hence most important. 
Hence, either method yields the same result. 

There are two special cases worth highlighting. The 
first concerns the 'scattering experiments' where the ini- 
tial data is a coherent compact pulse incident on the 
brane from the asymptotic region. We have ip(t,x) = 
f(t — x) before scattering, which allows us to swap dt for 
—d x in H67|) . We can also approximate ^^(x) « e lu ' x . 
Integrating by parts gives 



dxe iux ij;{Q,x). 



(68) 



Hence, X is the spatial Fourier transform of the initial 
data and \2\ 2 the initial power spectrum. We can use this 
to construct the Fourier transform of the brane signal, 
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keeping in mind that ipb(t < 0) = for this initial data: 



Z b (u') = I cfte-*" Vb(*) = J dte~ luJt ^ h (t) 

o 



it 6(tu — to ) + 



e l,5 T ? (w)J(w). 



To go from the first to second lines, we have used the 
Fourier transform of the Heaviside function. The integral 
associated with the first term in square brackets is trivial. 
The second integral expands as: 

oo 
1 (dx^x) j^ X ^\ 



The pole at u> = to' means that we must evaluate this 
as a principal value integral. Since x < 0, we complete 
the contour in the lower half plane where the scattering 
matrix is analytic (j) is analytic everywhere). Assuming 
the integral over the semi-circle at infinity vanishes, we 
obtain 



Z h {uj) = 2ire iS r ] (( J j)l(Lu), 



dE 



= 4^H|J( W )| 2 , (69) 



where we have defined the power spectrum in the usual 
way dE/dui — \Z b (oj)\ 2 . Physically, this says that the 
gravity wave power seen on the brane is given by the 
power in the initial pulse times £(w). In this sense, one 
could call £(w) a 'filter-factor', in that it tells us how the 
potential processes an input from infinity into a brane 
signal. Of course, influence of the QBSs is dominant 
in this process, because they give rise to the dominant 
features in £(w). 

The second special case concerns brane localized initial 
data: 



ip(0,x) = AS(x), ip(0,x) =BS(x), 
Re ir 



JM = 



2tt 



ioj 



(70) 



A similar argument to before gives the power spectrum 



dE 
duj 



A 1 



J 2 



(71) 



This is a little more complicated than the scattering ex- 
periment case (|69H in that knowledge of the trapping 
coefficient and initial data is not enough to determine 
dE/dw, one also needs to know R{uj). But we do know 
the behaviour near QBS resonances 



dE 
du) 



B 2 



A + — i 

LO z J 4(UJ — vo 



r 2 



(72) 



Hence, the system's resonances manifest themselves as 
Lorentzian peaks in the power spectrum. 



To summarize: In this section we have constructed 
the formal Green's function solution for the brane sig- 
nal. Both the scattering matrix and trapping coefficient 
appear explicitly in this solution, implying that the sig- 
nal will be dominated by QBS resonances. The brane 
power spectra arising from the scattering of compact 
pulses from infinity is completely described by the pulses' 
Fourier transform and £. However, more information 
is required for different types of initial data, including 
the brane-localized case. Regardless of the details, the 
prominence of rj within the Green's function ensures that 
dE/dto is peaked about the ES resonant frequencies for 
any type of initial data. 



B. Coherent states 

Knowledge of the power in a given signal is not enough 
to reconstruct the signal in the time domain; much in- 
formation is contained within the phase. So despite our 
understanding of the Green's function of the problem, it 
is still useful to conduct numeric scattering experiments 
to see what kind of behaviour is possible. 

In quantum potential problems featuring discrete 
wavef unctions, one can often construct 'coherent state' 
solutions that are wave packets whose dynamics mimic 
that of a classical particle travelling in the same poten- 
tial. In the current problem, we know that far brane 
configurations have a number of lightly damped QBS 
modes. For example, the (7,rt,,£) = (20,250,1) case 
has around seven modes with lifetimes ti/ 2 ^ 1000. On 
shorter timescales, these act like discrete bound states. 
Hence we can construct wave packets that behave like 
classical particles; for example, they can bounce back 
and forth within the resonant cavity. Of course, as time 
passes the QBS modes decay resulting in the decoherence 
of the pulse. 

To see these effects explicitly, we consider a scattering 
process in the ('y,r\ > ,L) — (20,250,1) case. The results 
of our numeric integration for 2p are shown in Fig. 1151 
Cosine modulated Gaussian initial data is fired at the 
potential barrier from well inside the photon sphere: a 
good percentage of this is reflected back to the black 
hole, while the transmitted part has had a redistribution 
of the frequencies which make it up, with preferential 
selection going to the frequencies corresponding to max- 
ima of the trapping coefficient rj. After scattering off 
the brane, the lowest frequencies in the signal become 
effectively trapped between the brane and the barrier, 
and bounce between them. While some leakage occurs 
at each bounce off the barrier, they cannot penetrate the 
barrier far, which means that they persist for a signifi- 
cant number of bounces. If we were to increase L in this 
experiment, we would see the bouncing last more or less 
indefinitely since the leakage at each bounce is so tiny. 

To a brane observer things can look pretty strange. If 
the initial pulse is particularly narrow, the bouncing of 
the signal between brane and barrier would result in a 
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FIG. 15: A cosine-modulated Gaussian pulse becomes tem- 
porarily trapped in the potential well for the case (7, rb, L) — 
(20,250,1). Initial data is chosen as = cos(a;) exp[— (x — 
xo) 2 /2a 2 ] centred at £0 = —20 with a variance of a = 5, 
moving to the right. In Fourier space this signal is a Gaus- 
sian centred about w = 1. We show the signal undergoing two 
reflections off the brane; the potential barrier at the photon- 
sphere partially reflects the signal resulting in a bouncing, 
temporary trapping of the wave. 



brane signal which disappears and then reappears many 
times over — a state which looks like a collection of recur- 
rent gravitons. This effect is shown in Fig. 1161 for a long 
time integration. 

The Fourier transform of this signal, Fig. 1171 shows 
which modes dominate the waveform. The main curve 
shown the FT of the brane signal; the Gaussian is the 
FT of the initial data, and the inset shows the signal 
at very late times. We can easily discern the peaks in 
the spectrum which correspond to the spikes in r], shown 
in the inset of Fig. Indeed, the signal in Fourier- 
space looks very like (initial data) x 77, as expected from 
our Green's function analysis. Note, however, the the 
lowest modes are far too narrow to resolve numerically 

(~ 10- 10 ). 
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FIG. 16: The brane signal for the situation shown in Fig. 1151 
The gravitons appear in short bursts in this model, despite 
arising from a single coherent state in the bulk. Only at late 
times does this morph into what looks like 'normal' behaviour 
for gravity waves on a brane. 




FIG. 17: Fourier transform of the brane signal in Fig. 1161 The 
original signal is a Gaussian in Fourier space (grey curve), 
which gets combed into discrete frequencies over successive 
oscillations between the photon sphere and the brane. The 
very late time part of the signal is shown in the inset: the 
higher frequency components have a quicker damping time so 
only the lowest QNMs are left. 



C. Steady state gravitational radiation: 
brane-black hole thermal equilibrium 

Our last 'initial value problem' example does not really 
involve initial data at all. Consider a situation where the 
gravity waves near the horizon are a super position of 
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standing waves: 



ip(t,x)= / du Zoaiu)^"* cos(ojx + 5), 



BOUND 
STATES 



(73) 



This represents a situation where there is an equal 
amount of energy flux from and to the black hole. Phys- 
ically, this is a good model of a brane in thermal equi- 
librium with the flux of Hawking radiation coming from 
the horizon. In that case, the power in the asymptotic 
waveform is given by the blackbody formula 



J/Tu - 1 1 



(74) 



where the Hawking temperature is Th = k/2tt and n 
is the surface gravity. Somewhat trivially, we use the 
definition of the trapping coefficient to obtain the brane 
signal 



and power spectrum 



dio r?(a;)Z 00 (a;)e l 



— = 4*«a,)|Z 00 ( W )| = ^ /Th _ 1 . 



(75) 



(76) 



This result is pretty much identical to what we had from 
the Green's function answer for the scattering of a com- 
pact pulse: the power of the brane signal is the power in 
the waveform at infinity times Hence, the same con- 

clusion holds here: The blackbody spectrum is amplified 
at the positions of QBS frequencies, and that amplifica- 
tion has a Lorentzian lineshape whose width is propor- 
tional to the lifetime of the resonance. 



VI. DISCUSSION 
A. Summary of resonances 

One of the main objectives of this paper was to gain 
complete knowledge of the quasi-normal resonances of 
Einstein-static braneworlds. Practically, this meant cal- 
culating QNM frequencies in different regions of the 3- 
dimensional (7, rb, L) parameter space characterizing the 
bulk gravity waves. Using two different numerical meth- 
ods and a variety of analytic techniques, we managed 
to determine the behaviour of the resonant modes in a 
number of different interesting limits. In Fig. 1181 we have 
summarized the key qualitative features of the resonances 
as a function of the model parameters. 

Another way of presenting our results involves classi- 
fying resonances as over-damped (Imw 3> Rew), moder- 
ately damped (Imw ~ Rew), or under-damped (Imw <C 
Re a;). Then we can ask: Which parameter choices pro- 
duce which type of resonance? 
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FIG. 18: Dominant resonant behaviour in different regions of 
the (7, Tb,L) parameter space 



Over-damped modes: Evanescent resonances occur 
when the brane is very close to the horizon - gravi- 
tons are sucked off the brane with great efficiency. 
In fact, when the brane is extremely close the con- 
cept of gravity waves is a bit dubious: all QNMs 
are purely imaginary, so no oscillations take place 
at all, and such modes just decay exponentially. 
For branes a bit further out there may be a small 
number QNMs with non-zero real part. 

Moderately damped modes: When the brane is lo- 
cated at a moderate distance from the black hole, 
or the black hole is large so that no resonant cav- 
ity exists, we have the situation of 'normal' QNMs. 
This means that any signal will typically be dom- 
inated by the mode with the lightest damping - 
which also has the lowest frequency. This mode 
tends to dominate any signal very quickly. As dis- 
cussed in detail in our previous work Q, all the 
resonant modes of pure tension branes are in this 
class. 

Under-damped modes: Very long-lived modes arise 
in several different situations: 

Zero-mode: For all (7, L), as n, — > 00 the 
fundamental mode becomes purely real, and 
matches the GR frequency; this mode is stuck 
to the brane. The higher overtones remain 
complex. 

Quasi-bound states: When the black hole is small 
(7 ^ 2) there exists a well in the potential; 
provided the brane is located past this well a 
pseudo-resonant cavity forms and QNMs with 
w 1 less than the peak in the potential may 
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be interpreted as quasi-bound states. These 
modes have tiny imaginary parts implying a 
very long lifetime before they tunnel into the 
horizon. 

Bound states: Quasi-bound states become bound 
on small scales: as L — > oo, we have Im(w) 
vanished exponentially quickly. This implies 
that on small scales an ES brane will look like 
a two brane model with virtually even spacing 
between all modes except the fundamental. 

This list and Fig. ^] encapsulate our essential results on 
the quasi- normal modes of Einstein-static braneworlds. 



B. The key role of the trapping coefficient in the 
initial value problem 

An intriguing, and somewhat unexpected, by-product 
of our work has been the trapping coefficient r\. Recall 
that this was originally introduced to help us calculate 
quasi-bound state frequencies when rb became too large 
for the series method to handle. But in our Green's func- 
tion analysis of the initial value problem in Sec. IV Al r\{io) 
played a prominent role Fourier transform of the brane 
signal. Indeed, for the case of a compact pulse propa- 
gating towards the brane from infinity, the power spec- 
trum on the brane was shown to be simply £ = rj 2 times 
the power spectrum of the pulse. In other words, the 
trapping coefficient tells us to what degree each Fourier 
component from x = — oo is localized on the brane. 

In this sense, the precise functional form of is cru- 
cial. To see how, assume for a moment that we only know 
the quasi-normal modes of the system and are ignorant of 
Then, in any given scattering experiment we would 
predict that the brane power spectrum will be enhanced 
in some range Aw ~ T n around each resonant frequency 
w n ; but we would not know beforehand to what degree 
each mode would be excited by the initial data. This 
is a known conundrum in black hole perturbation theory 
|3l| : quasi- normal modes by themselves only tell half the 
story, one also needs their 'excitation strengths' to really 
make detailed predictions. 

However, using £ to predict brane signal side-steps 
the whole issue because information about the excita- 
tion strengths is already encoded in the heights of the 
Lorentzian peaks representing the resonances. In small 
black hole/far brane cases £ is a very spiky function, 
which means that it effectively combs initial data into 
a sum of lightly damped sinusoids centered about quasi- 
bound state frequencies, and efficiently suppresses the 
power in-between. Even when the QNMs are heavily 
damped the trapping coefficient plays an important role, 
amplifying signals in the region near the fundamental 
QNM, but suppressing signals of lower frequency. At the 
other end of the spectrum, £ tells us in detail how high 
energy waves from the horizon are felt on the brane, a 



subject on which quasi-normal modes have nothing to 
say. 

But there are limitations to what can be learned in this 
fashion. If initial data has support outside the asymp- 
totic region (i.e., close to the brane), things become mud- 
dled because our simple Green's function result for the 
scattering of distant pulses does not apply. The only 
sure thing in these cases comes from the quasi-normal 
modes, in that they still predict the position and width 
of resonant peaks in the brane signal. The actual power 
carried in each mode must be obtained by other means. 
Having stated that caveat, the trapping coefficient re- 
mains a powerful tool for describing the resonances of the 
Einstein-static braneworld. It has the attractive quality 
of being much easier to calculate than QNMs from a se- 
ries solution, and it neatly communicates most of the im- 
portant features of the scattering problem. We expect it 
will be of considerable use in other 'one-sided' brane scat- 
tering problems that feature lightly damped resonances 
and a flat asymptotic region. 

C. Implications for cosmology and future work 

Einstein-static branes become closed cosmological 
branes if we allow their radii to depend on time rb = a(t). 
The dynamics of a(t) is governed by the Friedman equa- 
tion (|15fl . The evolution of gravity waves in the most 
general scenarios is more complicated that the static case 
considered in this paper, with a formal solution of the 
form 

ip h (t) = J dt' dx' G(t,x h ;t',x')X(t',x'). (77) 

Here, G is the real space retarded Green's function and 
X represents a generic source which includes ?/>(£, x) and 
perhaps other fields. Solving integro-differential equa- 
tions of this type are the essential obstacle in braneworld 
cosmology 0, |33- This is largely because the moving 
brane boundary condition makes it exceedingly difficult 
to deal with the Green's function in generic situations. 

However, in the high-frequency regime we can make 
use of a multiple-scales approximation |l7j | to make some 
progress. This relies on the fact that if gravitational 
waves oscillate much faster than the brane moves, they 
will 'see' the brane to be stationary. To zeroth order in 
this approximation, the Green's function above is just 
the frequency space Green's function we developed in 
Sec. transformed into the time domain. On dimen- 
sional grounds, 'slow' brane motion implies uj 3> H, 
which means that our static results should apply to waves 
with oscillation periods much shorter than the Hubble 
time. Another way to state this is that the static Green's 
function is a good approximation to the moving Green's 
function provided H(t — t 1 ) <C 1. 

Hence, any quasi-normal modes we have found will 
be approximate resonant solutions to the moving brane 
problem if Re uj n 3> H . In this picture, one thinks of the 
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frequency of such modes as quantities evolving in time in 
a quasi-static, adiabatic process. Then, plots like Figs. 01 
and 1101 which show the variation of uj n with fb, can be 
re-interpreted as depicting the evolution of resonant fre- 
quencies with cosmological epoch, provided we only pay 
attention to the high frequency modes. Ideally, we would 
like to understand the systems resonances at lower fre- 
quencies, which would involve moving beyond the zeroth 
order static approximation to the the Green's function. 
An interesting future project would involve formalizing 
such an approximation scheme and then defining gener- 
alized, dynamic quasi-normal modes. 

Finally, we note that due to the lack of time symmetry, 
it is more common in cosmology to characterize modes 
by their comoving wavelengths as opposed to their fre- 
quencies. For branes in the pure tension or far brane 
regimes, all resonances we found had Re w n > k, as might 
be expected from the null momentum condition in 5 di- 
mensions. Hence, our slow motion approximation can be 
re-cast in terms of spatial wavelengths: on sufficiently 
sub-Hubble scales, gravity waves on a cosmological brane 
will be dominated by the quasi-normal resonances of the 
corresponding static brane. Most importantly, this means 
that we recover a RS-like zero-mode for late cosmologi- 
cal epochs and on small scales. This is a critical test for 
the viability of braneworld cosmological modes involv- 
ing a bulk black hole - if it were not true, we would see 
significant departures from general relativity over astro- 
physical distances in the current universe. However, it is 
tantalizing that we only find an approximate zero mode. 
This leaves the door open to finding small corrections to 
standard gravity induced by the presence of a bulk black 
hole. This is an important subject for future study. 
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APPENDIX A: THE CLOSE BRANE LIMIT AND 
EVANESCENT MODES 

In this appendix, we find an analytic solution to the 
master wave equation when the brane is close to the black 
hole horizon r^ ~ 1. The goal is to obtain an exact ex- 
pression for the purely imaginary evanescent modes seen 
in Sec. lrTTCl for r h -> 1. 

The extreme close brane limit is defined by r^ — 1 <C 1. 
Since r G (l,rb), we have 



/«2«(r-l) 



e 2KX < 1, 



(Al) 



in between the brane and the black hole horizon (recall 
that k = (2 + j 2 )/"f 2 is the surface gravity). We expand 



the wave equation l|21(l to leading order in e 2KX , and ap- 
ply a linear coordinate transformation X = 2k (xb — x) 
to simplify matters. Assuming ipk = e luJt ®ku(X), and 
defining the parameters 



2k 2 \ r h 

2 _ [6 +QL 2 + 2L + 3)7 2 ](r b -l) 
2(2 + 7 2 ) 



(i 2 = 



(A2) 



yields the simple wave equation 

n 2 <f> ku (X) = -$L(X)+/? 2 e- x $ few (A), (A3a) 



$L(fj) = -A$ fc<1 ,(0). 



(A3b) 



Note that the X coordinate has been selected such that 
the brane is at X = and the horizon is at X — oo. The 
ODE i|A3a|) has the following exact solution satisfying 
~l3bl): 



®kUX) = F 1 I u (2Pe- x/2 ) + F 2 K^(2f3e- x/2 ), (A4) 
where v = 2iQ and 



Fx = pKx- v (2P) 
F 2 = (3h +u (2(3)l 



- (£0 + A)K v (2/3), (A5a) 
(iVL~ A)I V {2(3). (A5b) 



Here, J M and are the modified Bessel functions of 
order /j,. In order to find the quasi-normal frequencies, 
we need to know about the asymptotic behaviour of these 
Bessel functions near the horizon at X — oo. From the 
relations Hal 



I^z) = -r 



1 



(j, 



1 



+ 0(z 4 ) 



K^2z) = 



7r/_ M (2«)-J M (2z) 



2 sin7r/i 
we deduce the large X behaviour 



L(2pe 



-X/2} 



const, x e 



-mx 



(A6a) 
(A6b) 

(A7a) 



K- V {2f3e- X/2 ) -> const, x sin(f7A + tp), (A7b) 

where ip is some constant (real) phase. We see that the 
former represents a purely outgoing wave at the horizon, 
while the latter is a mixture of incoming and outgoing 
radiation. To have a QNM, we need to set the contribu- 
tion from the incoming radiation equal to zero. Hence, 
the quasi-normal frequencies are the (complex) to = 2kQ 
for which 



F 2 = F 2 {cu ;l ,r h ,L) = 0. 



(A8) 



Notice that r^ — 1 <C 1 means that (3 is a small pa- 
rameter. Furthermore, for any given choice of 7 and L, 
we see that A is of the same order as f} 2 . Therefore, it 
makes sense to expand F 2 in powers of (3. Again using 
I6all. we find 



F 2 =P 



1a 



/? 2 (i 



(1 

(1 + v)\ v\ 



(A9) 



22 



To leading order in (3, we find that F2 = implies Q/v\ = 
0, which has solutions 

= 0, v = -1,-2,-3... (A10) 

Expressing f2 and in terms of w and re yields: 

lim uj n = inn, n = 0, 1, 2, 3 . . . . (All) 

r b ->l 

This simple formula holds for all 7 and L. So, when the 
brane is very close to the horizon, the quasi-normal fre- 



quencies are evenly spaced along the positive imaginary 
axis, in agreement with the results obtained from the se- 
ries solution in Sec. IIII CI It is interesting to see that 
the QNM frequencies scale with the surface gravity of 
the black hole in this limit, just like the small-7 QNM 
frequencies of a S-AdS black hole |34|_or the overtones of 
a pure tension Einstein-static brane j3| • Finally, we have 
quantitatively compared (|A11|) to our series solutions for 
the cases shown in Fig. 0] and found excellent agreement. 
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